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Abstract. Wc study the effects of random scattcrcrs on the ground state of the one- 
dimensional Licb-Linigcr model of interacting bosons on the unit interval in the Gross- 
Pitaevskii regime. We prove that Bose Einstein condensation survives even a strong 
random potential with a high density of scatterers. The character of the wave func- 
tion of the condensate, however, depends in an essential way on the interplay between 
randomness and the strength of the two-body interaction. For low density of scattcrcrs 
or strong interactions the wave function extends over the whole interval. High density 
JZ* of scatterers and weak interaction, on the other hand, leads to localization of the wave 

function in a fragmented subset of the interval. 



1. Introduction 

While the effects of random potentials on single particle Schrodinger operators [I] and 
ideal Bose gases [2] [I] are rather well explored, the present understanding of such effects 
on many-body systems of interacting particles is much less complete. In recent years, 
however, many papers concerning the interplay of Bose-Einstein Condensation (BEC) 
and disorder have appeared of which references [S] [IS] are but a sample. 

In this paper we present results on a model that is in a sense the simplest one imaginable 
where this interplay can be studied by rigorous mathematical means. This is the one- 
dimensional Lieb-Liniger model [IT] of bosons with contact interaction in a 'flat' trap, 
augmented by an external random potential that is generated by Poisson distributed 
point scatterers of equal strength. We study the ground state and prove that, no matter 
the strength of the random potential, BEC is not destroyed by the random potential 
in the Gross-Pitaevskii (GP) limit where the particle number tends to infinity while the 
coupling parameter in a mean-field scaling stays fixed. The character of the wave function 
of the condensate, however, depends in an essential way on the relative size of the three 
parameters involved. These are the scaled coupling parameter 7 for the interaction among 
the particles, the density of the scatterers u, and the strength a of the scattering potential. 
All these parameters are assumed to be large in suitable units. 

The investigation consists of three parts. In the following Section 2 we prove Bose- 
Einstein condensation in the ground state of our model in a GP limit. This result holds, 
in fact, for quite general external potentials, provided they are bounded below. Specific 
properties of the random potential enter only in the estimate of the energy gap of an 
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effective mean-field Hamiltonian that sets a bound to the depletion of the condensate for 
finite particle number. 

Section 3 is concerned with the properties of the condensate in dependence of the 
parameters. In particular, for a ^> 1 we identify three different "phases" of the condensate. 
For large 7 ^> v 2 the condensate is extended over the whole trap (the unit interval in 
our model). A transition from a delocalized to a localized state takes place when 7 is of 
the order u 2 , in the sense that for 7 <C v 2 the density is essentially distributed among a 
fraction A 1 of the v 3> 1 intervals between the obstacles. For 7 3> v/{\av) 2 we still 
have \v ^> 1, but for 7 ~ u/ilnu) 2 the fraction of intervals that are significant occupied 
shrinks to 0(z/ _1 ). We stress, however, that in all cases there is complete BEC into a 
single state in the limit when the particle number tends to infinity. 

In an appendix, we prove an estimate on the energy gap of one-dimensional Schrodinger 
Hamiltonians. This estimate, applied to a mean field Hamiltonian, is needed to establish 
a lower bound for the condensate fraction that depends only on the parameters of the 
problem and not on the individual realizations of the random potential. The estimate 
on the energy gap is a generalization of a result from [20J and is of independent interest. 
Finally, for convenience of the reader, we collect in a second appendix some facts about 
the Poisson distribution of obstacles. 

Our study of the Gross-Pitaevskii energy functional in Section 3 can be regarded as 
complementary to the investigations in [5J, [7J El HOI [HI EH dS] where transitions between 
derealization and localization for solutions of the GP equation (also the time- dependent 
one) are considered from different points of view. Among these papers [7J is the one 
closest to ours as far as the questions asked are concerned, and it is appropriate to make 
a brief comparison. 

One difference is that the random potentials considered in [7J are intended to model 
experimentally generated laser speckles and somewhat different from the present model. 
In both cases, however, there is a length scale associated with the random potential (z/" 1 
in our case, <7r in [7J), and the strength of the random potential that is denoted by Vr 
in pj] is the analogue of our va . Moreover, our coupling parameter 7 may be compared 
with the chemical potential \i in [7J. 

The main difference, however, is that in [7J BEC is stated to hold only for large \i 
(analogous to large 7 in our model), while for weaker coupling a breakdown of BEC 
through a transition to a fragmented condensate and finally to a Lifschitz glass phase is 
claimed. We, however, 'prove rigorously that complete BEC holds in the whole parameter 
range considered when the particle number tends to infinity. Some further comments on 
this issue will be made at the end of Section 2. 

2. The Many-Body Model and BEC 

The model we consider is the Lieb-Liniger model of bosons with contact interaction on 
the unit interval, with an additional external potential V. In the next section we shall 
take V to be a sum of delta functions at random points in the unit interval, but for the 
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proof of the energy bounds and BEC in the present section the only assumption used is 
that V is a nonnegative sum of delta functions and a locally integrable function. The 
Hamiltonian on the Hilbert space L 2 ([0, 1], dz)®^™™ is 

TV 

H ■= E (~ d l + V (zi)) + ^ E *(* - ^) C 1 ) 

i=l i<j 

with 7 > 0, and we assume Dirichlet boundary conditions at the end points of the unit 
interval. The Hamiltonian has a unique ground state and we denote the ground state 
energy by E$ . 

The Gross-Pitaevskii energy functional is defined by 

£ GP ii'] := / 1 {|^(^)r + ^)l^)| 2 + (T/2)|^(^)| 4 }^. (2) 
Jo 

Again we assume Dirichlet boundary conditions. By standard methods one proves that 
there is a unique nonnegative minimizer ipo e H^([0, 1]), normalized so that J i/jq(z)(Iz = 1, 
with a corresponding energy e GP . This is equal to the ground state energy eo of the mean- 
field Hamiltonian 

h:=-d 2 z +V(z) +1 Mz) 2 -l J ^ (3) 
with ipo as the corresponding eigenf unction. 
THEOREM 2.1 (Energy bounds). 

E QM 

eo > > e (1 - (const.)iV- 1 / 3 min^ 1 / 2 , 7 }) . (4) 

Note that the error term depends only on 7 besides N and is independent of V (as long 
as V > 0). 

We denote the second eigenvalue of h by e\ and by N the average occupancy of ipo m 
the many-body ground state \l/o of H, i.e., 

iV :=trp«|^)^o| (5) 

where p 1 - 1 ^ is the one-particle reduced density matrix of ^o- The depletion of the conden- 
sate is (1 - 

THEOREM 2.2 (BEC). 

^1 _ ^ < ( const .)_^_iV-V3 min{7 l/2 )7} _ (6) 

Theorem 12.21 implies in particular complete BEC in the limit N — > 00 with 7 and V 
fixed. In fact, Eq. (jfJJ) implies that the reduced one-particle density matrix divided by N 
converges in trace norm to the projector on ipQ. The energy gap e± — eo is always strictly 
positive since h has a unique ground state. Its dependence on 7 and V is discussed in 
the Appendix. For the random external potentials introduced in (120]) below, Eq. ( 11031) 



4 



R. SEIRINGER, J. YNGVASON, AND V. A. ZAGREBNOV 



implies that complete BEC holds not only for almost every random sample of scatterers, 
but the depletion of the condensate goes to zero even in LP norm on the probability space 
for every p < oo. 

Our proof of BEC in the GP limit is simpler than the corresponding proof in three and 
two dimensions [18] because the one-dimensional case considered here corresponds to a 
high-density, mean-field limit. In contrast, the work in [18] deals with a dilute limit that 
requires quite different tools. 



2.1. Proof of Theorems 2.1 and 2.2. The upper bound in (TjJ is immediately obtained 
by using as a trial function for the expectation value of the Hamiltonian (pQ). 

For the proof of the lower bound in (j3J), as well as the proof of (jSJ), we begin by splitting 
off part of the kinetic energy and write, with e > and Pi := — id Zi , 

H = E { f 1 - ^r £ ) + n*) } + 4 E [fa + + - *)] • (?) 

i=l ^ ^ ' ' i<j 

Next we estimate the term in square brackets from below by a regular interaction potential 
using the bound 

- ed 2 z + 7<S(z) > w(z) (8) 

with 

w{z):= l + (l/2e)i eM - lzl/b) (9) 
for any b > 0, cf. [191 Lemma 6.3]. Note that 

5 b (z) := ^exp(-\z\/b) ^ 5(z) (10) 

as b — > 0, and J 6b(z)dz = 1 for all b. From (jSJ) we obtain the following inequality (in the 
sense of quadratic forms on wave functions that satisfy Dirichlet boundary conditions at 
the boundary of [0, 1]^ and are identically zero outside of [0, 1]^) 

N 

H * E { i 1 - 1) p" + + % E ~ *j) (") 

i=l i<j 

where we have denoted + (b^/2e)\ by 7 for short, and bounded — (N — l)e/2N from 
below by — e/2. Since 5b is of positive type we can, for any density p(z), bound the last 
term by a sum of one-body terms: 

2^ te -z,)> Wt^M^)-| f{S b *p)p-^z5 i (0)} (12) 

i<j i=l ^ ^ J 

where we have used that 7 > 7. We take p to be the GP density (and zero outside of 
[0, 1]). For the comparison of the mean field Hamiltonian ([3]) with the right side of (fT2"]) 
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we use the following simple estimates (recall that all integrals are effectively over the unit 
interval) : 

||^*p||oo<||^o||L<ll^o|| 2 ||^l|2<e; /2 (13) 

and 



\5 b * p- p\ 



oo = SUp 

X 



< 



sup 

x 




\p'(u)\duS b (y)dy 



x-y 



[p(x -y)- p(x)}5 b (y)dy 

<2\\MUeT J \y\ 1 '%{y)dy<2^^4! A b^. (14) 
Using 7 > 7 - &7 2 /(2e), Eqs. (p|)-(IIl|) and © lead to 



- \ 2 1 2e 2V2 ' 
i=i y 



ANb J 



(15) 



In the second term we estimate pf by e leading to an optimal e ~ 7e . Optimizing 
the resulting expression 



(const.)7 \b l/2 eT + 



over 6 gives b ~ N 2 ^ 3 e ^ 2 and 



E 



QM 



AT 



> e (l - (const. )A r " 1/3 7eo 1/2 ) . 



(16) 



(17) 



For large 7 we may use the bound e > (7/2) J p 2 > (7/2) while for 7 < 1 we use 
that eo > vr 2 since 7r 2 is the lowest eigenvalue of — d\ on the unit interval with Dirichlet 
boundary conditions. Altogether we obtain (J4]). 

To prove Theorem 12. 21 we trace the one-particle operator on the right side of (fT5|) with 
the one-particle reduced density matrix of the ground state of H and obtain the more 
precise lower bound 

N 



?Q M AT / AT \ 

> ^e + (l - ei - (const.) e m i n { 7 1 / 2 , 7} 



(18) 



where e\ is the energy of the first excited state of the mean field Hamiltonian h. Combined 
with the upper bound E^ M < Ncq this leads to the claimed estimate for the depletion 
of the condensate. 



2.2. Remarks. 1. If N<k denotes the occupation of the k lowest eigenvalues of h with 
energies up to e^, then the following generalization of ($5$) holds: 
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The proof is the same way as for ([6]), using again the bound ( fT5l) . For finite N the bounds 
dH]) and (TjUl) are, of course, only useful if their respective right-hand sides are less than 1. 
This may hold for (fl9|) even if it does not for ([6]). 

2. As remarked in the Introduction, the authors of [7J claim a transition from a BEC 
to a fragmented condensate and finally to a Lifschitz glass phase in their model. Since no 
rigorous limit of large particle numbers is considered in [7J there is no contradiction with 
the complete BEC that is proved in the present paper in the whole parameter range. In 
fact, for a given finite particle number N a fragmented condensate, i.e., (1 — N<k/N) small 
for some 1 <C k <C N, may be a reasonable substitute for the fully condensed state that 
emerges according to in the large TV limit. We also point out that as far as the density in 
the position variable is concerned, a fragmented condensate with non-overlapping single- 
particle wave functions is indistinguishable from a fully condensed state where the wave 
function of the condensate is a coherent superposition of the non-overlapping functions. 
The difference shows up, however, in the momentum distribution. 



3. The Gross-Pitaevskii Theory 

Having established the GP minimizer as the wave function of the condensate in the 
Gross-Pitaevskii limit of the many-body problem we now turn to the study of the depen- 
dence of this minimizer on a random external potential. Specifically, we shall take V in 
© to be 

V(z) = <tV u (z) (20) 

with a > and 

V u (z) :=^25(z-Zi/u) (21) 

i 

where v > and u is a variable in a probability space generating the random points Zi in 
R which are assumed to be Poisson distributed with density A = 1 (see Appendix 2). The 
distances ii between neighboring points Zijv and Zi+\/v are then independent random 
variables distributed according to the exponential law (cf. Eq. (11071) ) 

d Pu {£) = ue- £u d£. (22) 

Thus, with probability one only finitely many of these points are in [0, 1], and on average 
there are v such points. We start by considering the energy in an interval of length I 
between two of the points, which after a translation and scaling 

z -> x := z/£, a -> a := £a, 7 ->■ k := £7 (23) 

can be conveniently taken to be the unit interval. 

3.1. An Auxiliary Problem. For k > and a > 0, let e(n, a) denote the auxiliary GP 
energy 

e(«,a) := inf S K>a [</>] , (24) 
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where 

£*M ■= [ dx (l^'(^)l 2 + f l<K*)l 4 ) + f (I0(o)l 2 + I0(i)l 2 ) (25) 

for G i? 1 ([0, 1]). Standard methods show that there exists a minimizer for (124]) . which 
we denote by <f) K>a , i.e., 

e(«,a) = £ K , a [^K,a\ ■ (26) 

The minimizer is unique up to a constant phase factor, which can be chosen such that 
is non-negative. 

Note that, as an infimum over linear functions, e(n, a) is jointly concave in k and a. 
For us it will also be important that k H- Ke(/t,a) is strictly convex. This follows from 
the fact that 

















Jo 





«e(«,a)=inf <SiM\) ■ / M = K \ ( 27 ) 



together with the strict convexity of \<p\ 2 i— > £i >a 
We have the simple bounds 



and 



k f 1 K 

s(k, a) = S ,a(4> K ,a) + 2J dx I0«,a(a:)| 4 > e(0,a) + - (28) 

e(«,a) < £ K , a (0 o ,a) = e(0,a) + - y |0 o ,a(^)| 4 , (29) 

where we used the Schwarz inequality to obtain the first inequality. 

Obviously e(0, 0) = and e(0, oo) = n 2 . The minimizer 0o,a for k = is of the form 
(fio,a(x) = a a cos[b a (x — 1/2)], where b a > increases from to n as a increases from 
to oo, and a a > is determined by the normalization condition ||0o,q||2 = 1- A close 
inspection shows that ||0o,a|U is monotone increasing in a, and hence 

1 f 1 3 

dx |0 o ,aOr)| 4 < dx |0 o ,oo(a;)| 4 = - • (30) 
Jo * 

In particular, from (I28l) - fl30l) we obtain 

1 e(K,a) - e(0,a) 3 

2 " k " 4 [6l) 

for all k > and a > 0. 

Lemma 3.1. For some constant C independent of k and a we have 

e(«,oo) > e(/c,a) > e(«,oo) (l - Ca~ 1/2 ) . (32) 

Proof. The first inequality follows clearly from monotonicity of e(n, a) in a. To prove the 
second, it is enough to consider a > C 2 (for otherwise the right side of ( l32l) is negative). 
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Using a piecewise linear function which is constant for x G (e, 1 — e) as a trial function, 
we obtain (after optimizing over e) that 

e(«, cm) < | (1 + C«T 1/2 ) (33) 

for some C > and k large enough. In particular, since obviously e(n, a) > k/2, the 
second inequality in (1321 follows right away if k > a. 
We are left with the case a > k. Since 

e(/c,oo) > e(n,a) > ^ + a0 K , Q (O) 2 (34) 

(using the symmetry of <p K ^ a with respect to reflections at 1/2), we conclude from fl33|) 
that 

K , a (O) 2 < C— < CoT 1 ' 2 (35) 
a 

in this case. To obtain the desired bound, we use use A(0 K)O ,(x) — K , a (O)) as a trial function 
for e(n, oo), where A is an appropriate normalization factor. A simple calculation then 
yields @2). □ 

For later use we shall also introduce the Legendre transform of the map k h )■ Ke(n, a), 

g(n,a) := inf (ne(n,a) — yun) . (36) 

ra>0 

From the concavity of k H- e(«, a) and the uniqueness of the minimizer of (125]) one can 
easily conclude differentiability of the map n i— > e(n,a). Since k h- Ke(K,a) is strictly 
convex, the inffmum in (136|) is uniquely attained. If it is attained at some n > 0, this n 
satisfies the equation 

e(n, a) + n e'(n, a) = /x , (37) 

where we denote e'(/t, a) := d K e(n, a). We subtract the energy e(0,a) on both sides of 
(|37|) for convenience. Consequently, we observe that the infimum in (|36|) is attained at n 
satisfying 

n = [/i-e(0,a)] + — , : T^r^T • (38) 

e'(n, a) + -(e(n, a) — e(L), aj) 

Here [t] + stands for the positive part of a real number t. The unique solution of (138]) will 
be denoted by n(/i, a). 

From ( 13T]) and concavity of e(n, a) is follows that 1/2 < e'(K, a) < 3/4. Moreover, the 
second term in the denominator in ( 138]) lies in [1/2, 3/4], again by ( 13T]) . We thus conclude 
that the optimizing n satisfies 

2 

- \p - e(0, a)} + < n(/i, a) < \p - e(0, a)} + . (39) 

In particular, it is non-zero if and only if /i > e(0, a), and we can write 

n(/i, a) ~ [pi — e(0, a)], , (40) 
where a ~ b means that a/b is bounded from above and below by positive constants. 
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Finally, we note that since e(0, a) is increasing and concave in a, there exists a constant 
C > such that 

Ca 

eO,a>— -. (41) 
1 + a 

This bound will be useful later. 

3.1.1. The Average Energy. We shall consider the average energy for a particle distributed 
over intervals of side length £, whose distribution is given by the exponential law (I2"2"j) with 
v > 0. More precisely, we define 

e (7,z/)=inf f° dp„(£) ^-e(n(£)£j, oo) : i/ f°° dp„(£) n(£) = l\ , (42) 



where the infimum as over all n : M + — >■ K+ satisfying the normalization constraint. With 
the aid of f[3"oT) . we can alternatively write 

e (7, i/) = sup { fi + v I dp v {£) -^-g(fif, oo) I . (43) 
m>o I Jo £ 3 1 J 

From simple scaling it follows that 

eo(T^) =7eo(WV7)- (44) 

The infimum in ( 1421) is attained for n(£) = (£7) _1 n(/i^ 2 , oo) for suitable /i > 0, being equal 
to the optimal \i in ( j43|) . Using ( 1401) as well as the fact that 

1 < e x [ e~ l (t-—\dt<2 (45) 



we see that /x satisfies the relation 



t 



1 „ H e -*»/Vfi . (4 6 ) 
7 



In other words, \i ~ 7 f(i> 2 /'-f), where / : R+ — > M+ denotes the function 



1 for x < 1 . , 

(l+lnx) — 



Also e (7,f) ~ T/C^Vt)- 



3.2. The Gross-Pitaevskii Energy. We now turn to the actual GP energy functional 
we want to consider, 

CM := dz (V W| 2 + W*)| 2 + ||^)| 4 ) (48) 

with K, as in ( 12T1) . and define the GP energy as 

e w ( 7 ,<r,i/):= inf £ GP [0] . (49) 
<p 2=1 
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THEOREM 3.1 (Convergence of the energy). Assume that v — > oo, o — > oo and 7 — > 00 
in such a way that 

v v 

7 > 79 an d > \ — T-, ttt • (50) 

(lnz/) 2 1 + In (1 + z/ 2 / 7 ) 

Then, for almost every u, 

lim^#^ = l (51) 
e (7,^) 

where eo(7, v) is the deterministic function defined in fcfify , which satisfies eo(7, v) = 
je (l,v/^/) ~ 7/(^/7). 

3.2.1. Remark. Without the interaction, i.e., for 7 = 0, the energy is not deterministic 
but is simply the ground state energy in the random potential oV w . For o — > 00 it becomes 
equal to the kinetic energy of a particle localized in the largest subinterval free of obstacles 
and hence, with probability 1, of the order v 2 j (In v) 2 , cf. Appendix 2. 

We shall now investigate the optimal particle number distribution in (]42j) . Since 
n(fi£ 2 , 00) > if and only if fi£ 2 > ir 2 , the average number of intervals with non-zero 
occupation numbers is given by 

POO 

v / d Vv {t) = e - nu/ ^is=: X v . (52) 

Since v is the total number of available intervals, A < 1 in ( 1521) defines the fraction of 
them which are occupied. With the help of the definition of A in ([5"2]) the relation (146]) 
can also be written as 

71-21/2 T /rq\ 

= A* ~ T ■ (53) 



(In A" 1 )' A 
In other words, A is determined by 7 and v via the relation 

7 ~ ' (54) 

We can distinguish the following limiting cases: 

• If 7 ^> v 2 then by ( )54"|) we get A — > 1, i.e., all the intervals are occupied. The 
chemical potential satisfies /x ~ 7 in this regime. 

• If 7 ~ v 2 than A ~ 1, but A is strictly less than 1. Again we have \i ~ 7. 

• If 7 z/ 2 then A 1, i.e., only a small fraction of the intervals are occupied. The 
relation ( 153|) implies /i ~ (z//ln(z/ 2 /7)) 2 for the chemical potential. 

• If 7 ~ vj{\xYv) 2 then by (I54p the fraction A becomes 0(1/ v), i.e., only finitely 
many intervals are occupied. In this latter case, /1 ~ 7^ ~ z/ 2 /(lnz/) 2 , which 
corresponds exactly to the inverse of the square of the size of the largest interval, 
cf. Appendix 2. 
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In particular, A> 1/y only if 7 ^> z//(lnz/) 2 , and hence this condition guarantees that 
many intervals are occupied. In this case the law of large numbers applies and hence the 
energy becomes deterministic in the limit. If \v = 0(1), on the other hand, the value of 
e ui{li v) is random. This shows, in particular, that our condition (150]) on 7 is optimal. 

Also the second condition in (|50|) on a can be expected to be optimal. It can be 
rephrased as la ^> 1, where £ is the (weighted) average interval length 

POO 

l=v \ d Pl/ (£)£n(£) (55) 
Jo 

with n(£) = (£ / ~f)~ 1 n(fi£ 2 ,oo) the optimizer in ( H2l) . A simple calculation shows that 

l^U^il + XTlil + U 2 /^)). 

3.2.2. Remark. The conclusions above apply to the minimizing particle distribution in 
(H2j) . In the parameter regime defined in (150]) . they can be shown to apply also to the 
actual minimizer of the GP functional (EE8T) . This follows in a standard way by perturbing 
the energy functional, and then taking a derivative with respect to the perturbation. 
Specifically, one can add a term (3 Lf(\i[)(x)\)dx for suitable functions / to the GP 
functional, and change the definition of eo in (I42p accordingly. One then concludes that 
the energy asymptotics (1511) continues to hold in this perturbed case. Taking a derivative 
with respect to (3 at (3 = yields the desired information on the GP minimizer. 



3.3. Proof of Theorem EUJ 

3.3.1. Preliminaries. Let z\ < Z2 < • • • < z m denote those points in (121!) which lie in the 
open interval (0,z/). Let also z = and z m+ i = u, and denote £{ = \z i+ i — Zi\/v for 
< i < m. Then Y^=o^j = 1- For any ip G -f^o([0, 1]), we have 

m 

^ GP M = ES^^<,[^] (56) 
3=0 c o 

with 

rij = + dx\ip(x)\ 2 (57) 

and 

^(x) = \ 0j + tjx) for x G [0, 1]. (58) 

We can choose if} G Hq([0, 1]) such that ipj oc <^ 7;00 f° r & H j an d { n j} is an arbitrary 
collection of non- negative numbers adding up to 1. In particular, we obtain the upper 
bound 

m 

eu,(7, o-, v) < inf ^ if e (^7, 00) , (59) 
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where the infimum is over all collections of non-negative numbers {no,ni, . . . ,n m } with 
From (|56|) we also have the lower bound 

m 

eu,(l, <r, v) > inf -Aefatjl, V) ■ ( 60 ) 
|nj} j=o l i 

From Lemma \3. II we conclude that the two bounds (ISP]) and (1591) agree, to leading order, 
as long as IjO is large for all the intervals of length £j. In particular, the study of the GP 
energy reduces to a study of independent intervals in this case. More precisely, it suffices 
that tjO be large for those intervals containing most of the particles. 

3.3.2. Poisson Distribution. In this subsection we shall collect a few useful facts about the 
Poisson distribution, we will be used later on. A more detailed discussion of some of the 
points can be found in Appendix 2. Let Xm{w) denote the characteristic function of the 
set corresponding to having exactly m scatterers inside the interval (0, v). Its expectation 
value equals 

<*»> = ^ ■ (61) 

We shall frequently use the bounds 

e -" Y" — < e ~ v V ^-X m - Xu = e — (i-A+AinA) for A > ! ( 62 ) 

and 



ml A — ' ml 

m>\v m>0 



e^Y^-< e - u Y ——\ m ~ Xu = e -(i-A+AinA) for x<1 (63) 
^— ' ml ^— ' ml 

m<Xu m>0 

Note that 1 — A + AlnA>0forA>0, with equality only for A = 1. 

For any function depending only on £j for some fixed 1 < j < m — 1, we have 

POO 

(/&)) = / dp v (£)f(£) (64) 
Jo 

with dp v defined in ( 122|) . This fact (cf. Eq. (j!07p in Appendix 2) will be used repeat- 
edly below. Moreover, the variables £j and £k for k ^ j are independent, and hence 
(/i(^)/ 2 (4)> = (/i(^))(/ 2 (4)) in this case. 

3.3.3. Upper Bound. Let fi = ^(7, be the optimizer in (143]) . With n(/x, a) defined 
above, we then have 

/•oo 1 

" / t-™(^> 00) = 1 . (65) 



In every interval of length £j, 1 < j < m — 1, we shall place n.,- = (^•7) _1 n( / u^, 00) 
particles. For simplicity, we shall not place any particles in the first and last interval. The 
total number of particles in the system is then N = Yl^i n j- 



DISORDERED BOSE EINSTEIN CONDENSATES 13 

Lemma 3.2. If v — > oo and 7—7-00 with 7 ^> u/(lnu) 2 then 

limA^=l (66) 

for almost every u. 

Proof. We shall show that (N) — > 1 and (iV 2 ) — > 1 in the limit under consideration, which 
implies the result. 

With Xm(w) denoting the characteristic function of the set corresponding to having 
exactly m scatterers inside the interval (0, u), we have, by symmetry, 

(iV) = - VVm - l)< X ^r^(K, 00)) . (67) 

ry r 

' m>2 

To obtain a lower bound, we can restrict the sum torn > 1 + (1 — e)v, for some < e < 1. 
We thus obtain 

' m>l+(l-e> 

= (1 - e) 1 - - Yl (Xmi^n^il 00)) , (68) 

y ^ m<l+(l-e)v / 

where we used (J65l) the obtain the last equality. From ( J39l) it follows that n(/x£ 2 , 00) < /i£ 2 , 
and hence 

(JV)>(l-e)(l-^w) , P{»)= E XmH- (69) 



m<l+(l-e)f 



A Schwarz inequality yields 



(P4> < (P) 1/2 (f 2 > 1/2 = (P> 1/2 ^ , (70) 

where we have evaluated the expectation value of l\ using (IMj) . Moreover, it follows from 
(163]) that (P) is bounded by a factor that decays exponentially in v. In particular, using 
flU, it follows that 

liminf(A0 > 1 -e. (71) 
Since e was arbitrary, this proves the desired lower bound. 

For an upper bound, we can proceed similarly. From ( 1671) . ( |65l) and n(/i£ 2 , oo) < /jlf 
we have 

(N)<l + e+^ J2 (m-ViXmti) (72) 

^ m>l+(l+£)^ 

for e > 0. An similar analysis as above, using the Schwarz inequality as in ([TP]) , and (|62|) 
instead of (|63|) . then yields 

limsup(JV) < 1 + e. (73) 
This shows that (N) — > 1, as claimed. 
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Next we consider (N 2 ). Similarly to (l67j) . we can write it as 

W = ^E^- 1 )^ -^(XmtlM^l 00)4^(^,00)) 



^ m>3 



+ A E( m - l)(x m ^ 2 ^(^?, oo) 2 } . (74) 

^ m>2 

An analysis as above shows that the first sum can be restricted to an interval m G 
[(1 — e)v, (1 + s)v\, and therefore yields 1 in the limit considered. Hence it remains to 
show that last term in (1741) goes to zero. If we restrict the sum to m — 1 < 2u, say, then 
we have 

(m - l)( Xm ii 2 n{fiei, oo) 2 ) < 2u(£^ 2 n{fi£l, oo) 2 ) < i^e'™ 1 ^ (75) 

m=2 

where we used the upper bound in (139]) as well as the bound 



e x J e - * ft - y J dt < fc! 2 fc (76) 

for k = 2. It follows from (j4"6"]) that (1751) . when multiplied by 7~ 2 , goes to zero as z/ — >• 00 
if A > z//(lnz/) 2 . 

For the remaining terms corresponding to m > 1 + 2u, we use again the Schwarz 
inequality, which implies that 

(m- l)(x m ^«(K,oo) 2 ) < uV2e^- hl ^(£^n(^loo)Y 2 - (77) 

m>l+2u 

Because of the upper bound in (1391) and (176]) for A; = 4 the latter expectation value is 
above bounded by 

(£^n(^l oo) 4 ) < 4! 2 4 4 e-WVF . (73) 

With the aid of ( 1461) one now readily checks that the right side of ( 1771) . when divided by 
7 2 , goes to zero as v — > 00. This completes the proof. □ 

We can take rij = (£j7) _1 n(/i£ 2 , 00) for 1 < j < m — 1 as above, and divide each rij by 
N = l^Jj=i Tij in order to have the right particle number. Then ( 159]) implies the upper 
bound 

m— 1 

e„( 7 , a,u) < E = ^ ^ e ( n / iT /iV, 00) . (79) 

The following lemma concludes the proof of the upper bound in Theorem 13.11 
Lemma 3.3. If v — > 00 and 7—7-00 with 7 ^> u/(lnu) 2 then 

E , s 

lim — = 1 (80) 



00 



„2\ fc 
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for almost every u. 

Proof. We proceed as in the proof of Lemma [3.21 and show that both (E)/eo( , y,v) and 
(E 2 ) j eo(7 , v) 2 converge to 1 in the limit considered. Because of the result in Lemma \3. 21 
we can replace N by 1 without loss of generality. Using (|40l) (for a = oo) and (l3Tj) we see 
that 

-|e(n/j7, oo) ~ firij . (81) 

Since also /i ~ eo(7, v) we can proceed exactly as in the proof of Lemma f3T2l to conclude 
the result. □ 

We thus conclude that, for almost every u, 

,. e w (7,a,z/) 

hmsup — < 1 (82) 

eo(7>^) 

in the limit ^ —> oo, 7 — )■ 00 with 7 ^> ^/(lnz/) 2 . 

3.3.4. Lower Bound. We start from (1601) . Recall the definition of g(/i, a) in (1361) . We have 

m m - 

e u {l, 0-, v) > inf -f e(n.,-f,-7, ^cr) > /i + ^ —^g(fi£ 2 , Ip) (83) 
lTy} i=o i=o 7 *i 

for any pel. 

Using the lower bound in (}3~Ij) . as well as ( HIT) , have 

<?(^ 2 » > -1 [/i£ 2 -e(0,^ > -^-6(£-£), (84) 



where 



5) 



7 1 10 
" = -Ya + ^ 2 + ~, 

and # denotes the Heaviside step function. For the first and last interval, corresponding 
to j = and j = m, respectively, we shall simply use that 



7 ^K>V)>-^- (86) 



J 

When divided by \i ~ eg (7, f), this goes to zero almost surely in the limit considered, 
since £j < v~ x and ji <C uj. 

Pick some small e > 0, and consider the contribution to the sum in flHBl coming from 
intervals with length £j < (ea) -1 , for 1 < j < m — 1. If eo~ > l/£, this contribution is 
zero, hence we can assume ea < \ jl from now on, which is equivalent to 



2 „ + 



^<^- (87) 
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Because of our assumption ( 1501) on a, this means that 

v 



1 + In (1 + 1/2/7) 
i.e., v 2 <C 7. We have 

\g{ l il\Lo)> --^ ^— — ^ ^7— - (89) 

7 f 2 7 ea ^o-(l + ln(l + z/ 2 /7)) v y 

in this case, and the last fraction goes to zero in the limit considered. 

Finally consider the intervals 1 < j ' < m — 1 with length £j > (sa)^ 1 . In this case, we 
can use Lemma 13.11 to bound 

SK 2 ,V) > {l-Ce^g^l-Ce^r 1 ,^). (90) 

An either case, we thus have the lower bound 

1 (l - Ce 1/2 ) C11 

^ 9 (K, V) > ^r-kWb - Ce"r',oo) - ff(1+ ,„ (1 % /7)) PI) 

for 1 < j < m — 1, for some C > 0. 

Let /i = /i(l — Ce 1 / 2 )" 1 . With n(/i, 00) defined in Subsection 13.11 as the optimizer in 
(1361). we can write 



g(fx£ 2 , 00) = n(fi£ 2 , 00) e(n(fi£ 2 , 00), 00) — [i£ 2 n(fx£ 2 , 00) . (92) 
We shall choose /Z as in Subsection 13.3.31 i.e., such that 



f 00 1 

/ dp v {i)— n(fi£ 2 , 00) = 1. (93) 
Jo «7 



Then /2 ~ e (-y,u). Let ^ = [£^y x n{ji£ 2 , 00) and iV = Y0j=\ n r From ®> flE 
and (1221) we have 



m— 1 



e w (7, ^ f) > (1 - C^ 1/2 ) U (1 - iV) + ^ -£e( nj £ jlt 00) 

C/im /j, 2 (£ + £ m ) 



(94) 



a(l + ln(l + 1/V7)) 2 7 

We have already shown in Lemma 13.21 that lim N = 1 almost surely, and in Lemma 13.31 
that the ratio of the sum in (f9~4l) and eo(7, v) converges to 1 almost surely. Moreover, 
after division by eo(7, v) ~ \i the terms on the last line of (1941) vanish almost surely. This 
proves that 

liminf^>l-C^ (95) 
for almost every u. This holds for all e > 0, hence the proof of the theorem is complete. 
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4. Conclusions 

The main conclusions that can be drawn from the preceding discussion are as follows: 

1. BEC in the ground state of the interacting gas in the Gross-Pitaevskii regime can 
survive even in a strong random potential with a high density of scatterers. As far as 
BEC is concerned the interacting gas in this regime thus behaves in a similar way as an 
ideal gas at zero temperature. The character of the wave function of the condensate, 
however, is strongly affected by the interaction. 

2. A random potential may lead to localization of the wave function of the condensate, 
even though the density of obstacles is much less than the particle density. The inter- 
particle interaction counteracts this effect and can lead to complete derealization if the 
interaction is strong enough. 

3. In terms of the interaction strength, 7, and density of scatterers, u, the transition 
between localization and derealization occurs in the model considered when 7 ~ v 2 . In 
the localized regime 7 <C v 2 the chemical potential is fi ~ [yj ln(V 2 / lr ))) 2 . For 7 < *V(ln v) 
the condensate is localized in a small number of random intervals. 

5. Appendix 1: Energy Gaps 

We consider the Schrodinger operator H = —d 2 + W(z) on L 2 ([0, 1]), with Dirichlet 
boundary conditions. Assume, for the moment, that W is a bounded and continuous 
function. Without loss of generality, we may assume W > 0. Let < e$ < e\ denote the 
lowest two eigenvalues of H . 

Lemma 5.1. Define rj > by 

rf = 7T 2 + 3 f W{z)dz (96) 
Jo 

Then 

ei-e > r/ln (l + 7re~ 2 ' 7 ) (97) 

For the proof, we follow closely |20j . 

Proof. For general E > 0, let u(z, E) denote the solution of 

- u" + Wu = Eu (98) 

with u(0, E) = and u'(0, E) = 1 (the prime denoting d z ). Pick an 77 > 0, and introduce 
the Priifer variables r(z, E) and 6(z, E) via 

u = r cos 9 , u' = —r]r sin 9 

The variable 9 is only determined modulo 2n. We fix it uniquely be requiring it to be 
continuous and 9(0, E) = —tt/2. We note that 9 is increasing in E for fixed z, i.e., 
<f(z, E) = d E 9(z, E) > 0. This follows from 
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where the dot stands for 8e- The function x — uu' — uv! satisfies x' — u2 an d is thus 
positive since it vanishes at z — 0. 
From dSBD it follows that 

0' = v ~ l (E- W) cos 2 6 + rj sin 2 6 (99) 

In particular, 6' > near points z G [0, 1] where u vanishes. Since u(z, eo) has no zeroes 
for z G (0, 1), we conclude that 9(1, eo) = tt/2. Moreover, u(z, e\) has exactly one zero in 
(0, 1), hence 6(1, ei) = 3ir/2. In particular, 

9(l,e 1 )-e(l,e )=* (100) 

We have </?(0, E) = for all £. From (JHSJ) it follows that 

<p'=(ri + Tf 1 (W - E)) sm(26)<p + r/" 1 cos 2 (0) 

< | + r]~ l (W - E)\ + 77" 1 

From Gronwall's inequality [211 Thm. III. 1.1], we thus conclude that 

<p(z,E)<~expU Q \v + V~ 1 (W(y)-E)\dy^j 

In particular, 

<p(l, E) < i exp ^77 + t/~^ E + 77" 1 jf 

for £ > 0. 

This implies that 

7r = 0(l,ei)-0(l,eo)= ! \(l,E)dE 

< exp (^77 + 77- 1 J W + 77-^0^ (exp (77- 1 (e 1 - e )) - l) (101) 

Using the trial function y/2 sin(7rz), we see that e < n 2 + 2 W. From f llOip we thus 
have 



ei — e > 77 In 1 + 7r exp —77 



^ fw(y)dy-^-\] 
V Jo V ) ) 



We choose 77 in order to make the exponent as large as possible, which leads to the choice 
(EEL and concludes the proof. □ 



By a simple approximation argument, the bound (|9"T|) extends to our case of interest, 
where 

m 

W(z) = aJ2$(z- Zi) + 7^o0?)| 2 (102) 

i=l 

for er > and 7 > 0. The result is then 



ei - e > 77 In (l + e~ 2,? ) with 77 = ^/n 2 + 3mcx + 37 . (103) 
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5.0. 5. Remark. While the bound ( 1 1 3 f> on the gap is uniform in the location of the scat- 
terers, we do not expect it to be optimal for large n. In fact, the gap is presumably 
smallest when all the equal to 1/2. In this simple calculation (for 7 = 0) 
shows that the gap is asymptotically equal to 32-7T 2 / (am) for large am. 

6. Appendix 2: Remarks on the Poisson Random Point Field (RPF) on the 

Real Axis R. 

6.1. Remark 1. (Poisson and exponential distributions, see, e.g., [3].) 

The following theorem is a consequence of the hypothesis of an independent uniform 
distribution of m point "impurities" on the interval A/ = [—1/2,1/2] C K (homogeneous 
m-point binomial RPF) and the limit: Aj — > R, m — > oo, for a fixed density 

A = lim — . (104) 

m, I— >oo I 

Theorem 6.1. (a) In the limit (j!04p the finite-volume binomial RPF of impurities {Xf} 
converges (in distribution) to the Poisson RPF {X u } with the Poisson distribution 

P{w : iXTlAJ =n} = e~ Xa , (105) 

of the number of impurities in any interval A a of the length a. 

(b) The uniform distribution of m — 1 independent points of impurities split up the 
box Ai into a set of sub-intervals -_ x with U™ =1 /" = A; and their random lengths 
{|/^| = L^YjLn u e ^ have the joint probability distribution 

dP l>m (L x , L m ) = V — x J ' 5(Lt + ... + L m -l) dL x dL 2 ...dL m . (106) 
Therefore, the random variables are dependent. 

(c) In the thermodynamic limit: A = lim^z^oo m/l, the intervals form an 
infinite set of independent random variables and the distribution corresponding to (I106|) 
converges (weakly) to the product-measure distribution a\ defined by the set of consistent 
marginals: 

k 

da XJl ,... tjk (L jl} .... L jf ) X k J] < x '-~dL h ,*6N. (107) 

s=l 



6.2. Remark 2. (Log intervals [3] , [3] .) 

We can estimate the maximal length, (maxi^^ m L")i, of an interval in the family {/j^} 
of random sub-intervals of A; = [—1/2, 1/2} when I — > oo. To this end we define the events 



_ (maxi<; Km L^)i ^ 5- 



A, := {co : 7, J/ < -} , / G N and I > 2 , (108) 

la I A 



20 R. SEIRINGER, J. YNGVASON, AND V. A. ZAGREBNOV 

as well as the following "tail" event: 

oo oo 

^:=limsupAf=f|U A f- ( 109 ) 

I— >oo -, , 

n=l l=n 

By definition, (I109P is the event that for increasing I the length of the largest subinterval 
in A; is bounded from above by the nonrandom number (<51nZ)/A for all but finitely many 
I G N. The following theorem holds: 

Theorem 6.2. Let A > be a mean concentration of the point Poisson "impurities" on 
R. Then 

(i) for any S > 4, one has: 

P«) = 1 , (110) 

(ii) for any 5 < 2, one has: 

P«) = 0. (Ill) 
Proof: (i) The complementary event to fllOSp is 

(AfY := {a, : ( max > (in/}. 

Let iVj := [A//(<51nZ)] + 1, where [x] denotes the integer part of x > 0, and define a new 
interval: 

Ai:=[-JVi(^lnZ), JV,(^ln01 => A i ■ ( 112 ) 

Split this bigger interval into 2Ni identical disjoint intervals {</'}^ = i of size <5(2A) _1 hiZ. If 
the event (Af) c occurs, then there exists at least one empty interval J\ (interval without 
any impurities), i.e., 

{A\) c C |J {u : J\ is empty} . (113) 

Note that for a given number of m — 1 impurities on the finite interval A; the RPF {Xf} 
is binomial for any two sets: J\ and A; \ J\. Since the probability for the interval J\ to be 
empty depends only on its size £(2A) _1 mZ and the size of A;, one obtains by (11121) and 
(11131) the estimate 

^-iM 5 i) C ) <m[l - — lnZ/|A,| < 2N t 1 - — . (114) 



2A " 7 V 2N i, 

To get the last inequality we put m — 1 = [A/] + 1 > A/, see (1104[) . and we use that 
{Ni - 1) Skil < XL Since the estimate ffTl4l) yields 

lim P^+liCW) 6 )/ [2([A//(51n/)] + 1) Z~ 5/2 ] < 1 , (115) 

we obtain that for 5 > 4 

E p iaW(^) c ) < 00 • ( 116 ) 
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By Theorem 16.11 binomial distribution Pmn +1 ;(") converges to the Poisson distribution 
P(-) on R, when I — > oo. Then estimate (IllSp and ( I116P imply 



< oo. 
1>1 



Therefore, by the Borel-Cantelli lemma P((y4^ c ) c ) = 0, which is equivalent to ( IllOp . This 
also implies that there exists a subset Q C f2 of full measure, P(f2) = 1, such that for each 
uj G Q one can find lo(u) < oo with 

r 

P {uo : ( max L?), ^ - In/} = 1 , for 5 > 4 , 

T<Km J A 

for all / ^ /o( w ) an d m > [A/]. 

(ii) Consider the event (Af ) c and note that one can estimate its binomial probability from 
below taking only one term in the union (11 13)) . Then 

^-m((^) c ) > U-^lnJ/MJ > (l--ln//[A/]J , (117) 

where we put m — 1 — [XI], see (11041) , and we use (11121) . Then ( 11171) yields 

hmP^((Af) c )/(r 5 / 2 )>l, (118) 

By the same arguments as in (i) one concludes from ( 11181) that for 5 < 2 



lf) c ) = oo. 



Since the events {(Af) c }; >2 are independent, by the Borel-Cantelli lemma we obtain: 
P((A^) C ) = 1, or equivalents (jTTTl) . This again can be translated as 

P {w : ( max L?), > - In/} = 1 , for 5 < 2 , 

for all / ^ /q(w) and m > [XI]- □ 



Since the Poisson RPF on K conditional to the interval [—1/2, 1/2] and to m — 1 impu- 
rities coincides with the binomial RPF, instead of increasing boxes with finite binomial 
RPF configurations one can consider restrictions of infinite Poisson RPF configurations 
to the increasing set of windows [— 1/2,1/2], cf [I]. Then Theorem 16.21 claims that almost 
surely (with respect to the Poisson RPF distribution) the maximal length of impurity- 
free sub-intervals, (max 1? ;j^ m ) h in this windows is bounded from below and above by 
(const) In //A. 
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